Antiphase synchronization of two nonidentical pendulums 
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We numerically study the synchronization of two nonidentical pendulum motions, pivoting on a 
common movable frame in the point of view of the dynamic phase transition. When the difference 
in the pendulum lengths is not too large, it is shown that the system settles down into the dynamic 
state of the antiphase synchronization with the phase difference vr. We observe that there is a 
bistable region where either the antiphase synchronized state or the desynchronized state can be 
stabilized. We also find that there exists a hysteresis effect around the dynamic phase transition as 
the length difference is adiabatically changed. 
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The phrase of odd kind of sympathy (in short, odd sym- 
pathy) was used by a prominent Dutch mathematician 
and physicist Christiaan Huygens in order to mention an 
interesting phenomenon observed for the two pendulum 
clocks attached on the wall [l| : Even after an intentional 
disturbance, two pendulum clocks evolve into a state in 
which they swing together in the opposite directions in 
synchrony. This is one of the most historical observa- 
tions which naturally lead to the concept of synchro- 
nization of dynamic variables, and have attracted many 
scientific/engineering researchers for a long time [l|, y]- 
In the language of synchronization study, the odd sym- 
pathy corresponds to the antiphase synchronization in 
which the frequencies of the pendulums' oscillations be- 
come identical but the phases show the mismatch tt. 

There are abundant examples of synchronization phe- 
nomena from biological objects ^ and celestial sys- 
tems [J] in nature to manmade electrical and mechanical 
systems @- Interestingly, such synchronization behav- 
iors have been in most cases empirically discovered by 
chance a posteriori (see the Introduction of Rcf . 2) , and 
we still need to understand more to predict why and when 
it happens. The lack of the precise knowledge can result 
in a noncontroUable outcome or a catastrophic disaster 
in the worst case Q . Even for such a simple setup of two 
pendulums, which allowed a chance for C. Huygens to be- 
come aware of the odd sympathy, investigations have still 
been performed to unveil the synchronization property 
in it 0, ITHIII . Such a study on the prototypical setup is 
necessary for a comprehensive understanding of the syn- 
chronization phenomenon as a dynamic phase. This is 
also to contribute to making the abstract-model-based 
synchronization research more fruitful in the viewpoint 
of the synchronization as a thermodynamic phase [13]. 

More than two hundred years later after the Huygens' 
observation of the odd sympathy, a qualitative under- 



standing was firstly tried by Korteweg in 1906 [3|- Nearly 
one hundred more years later, in early 2000s, Bennett 
et. al. revisited this issue [7'] and successfully repro- 
duced the antiphase synchronization in their experiment 
and theoretical model study. Other researchers of the 
two-pendulum system have also been interested in the 
inphase synchronization, in which the frequencies of the 
pendulums' oscillations are identical without the phase 
mismatch. The possibility of the existence of an inphase 
state was suggested in Ref. y, and was credibly repro- 
duced for the first time in the experiment introduced in 
Ref. S These works have been followed by a series of 
studies where various experimental setups and mechani- 
cal models have been proposed to understand the onset of 
inphase synchronization as well as antiphase one [9l4llj. 
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The motivation of this work is to investigate how 
generic the synchronization is in the two pendulum sys- 
tem. This is attributed to the fact that each work 
above jTHllj is based on its own specific experimental 
setup, for example, escapement mechanism and system 
structure. We are also interested in the bifurcation prop- 
erty of the synchronization state. In this paper, we nu- 
merically study the antiphase synchronization of the two 
nonidentical pendulums in the viewpoint of the dynamic 
phase transition. For the curiosity of generic feature, we 
introduce a model whose detail is as simple as possible. 
As the lengths of pendulums are varied, it is observed 
that there exists a broad range of bistable region, where 
the antiphase synchronized state and the desynchronized 
state coexist depending on initial conditions. The bista- 
bility is also shown to lead to a hysteresis effect as the 
lengths of pendulums are varied adiabatically. 



We first introduce the three degrees of freedom 
model lHI, schematically shown in Fig. [TJ We apply the 
Lagrangian least action principle ^15i] for nonconservative 
forces, i.e., the escapement and the friction forces in this 
work, and achieve the equations of motion in dimension- 
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FIG. 1: Each pendulum (i = 1 and 2) of the mass m is 
connected to the rigid common frame of the mass M by the 
massless rigid rod of the length li. The frame at the horizontal 
position X is attached to a spring with the spring constant 
K. It is assumed that the sliding motion of the frame and 
the pivot motion of the pendulum is dissipative, which are 
described by the friction constant B and 6, respectively. The 
hatched region stands for the rigid immobile wall considered 
as the reference frame and the uniform gravitational field g is 
applied vertically. The three movable objects (the frame and 
the two pendulums) are described by the coordinates X, 9i, 
and 02, respectively, which are the three degrees of freedom 
in the system. 



less form (see FigUand compare with Ref. 7) 

e, + 2j0, + (sin 9, + x cos 0,)/k - .h = 0, (1) 
i + 2ri; + n^x + mXI ^~»(^* <=°^^» ^ ^^ ^^^^*) = 0' (2) 



where the frame coordinate x{= X/l) and the dimension- 
less length li(= U/l) are measured in units of ^ [li — l+Ci, 
where e^ can be interpreted as a small but unavoidable 
relative error in manufacturing or measurement), and the 
time t (the dots on the symbols represent time deriva- 
tives) is in units of \jTfg, respectively. We have also 
defined the reduced mass as /i = m/{M + 2m), the effec- 
tive coupling strength of the frame as Vt^ = K/{M + 2m), 
and the dimcnsionless friction coefficients for the frame 
r = {B/2){yTfg)/{M + 2m) and for the pendulum 
7 = {b/2)^Jl/g, respectively. Since every degree of free- 
dom in the system is subject to the damping, the motion 
will eventually stop at x = ^i = 02 = in the absence of 
the external energy source. As the escapement method, 
we apply the impulsive force /, in Eq. ([Ij when 0^ = 
[fi < for 61^ > and fi>0 for 6i < 0] . Our escapement 
method enforces the pendulums not to stop by inject- 
ing kinetic energy into the system. It should be noted 
that the death phase observed in Ref. Ill is not allowed in 
our setup. Although other escapement mechanisms could 
lead to different results, we believe that our escapement 
algorithm is not too far from the reality, and that the 
generic qualitative results obtained in this work should 
be observable in properly prepared real experiments. We 
also remark that the other mechanical detail such as fric- 
tion in joints brings about a complex phenomenon, which 
has been thoroughly studied in Ref. |12| ) 

In the numerical experiments, we use 7 = 3 x 10^^, fl — 
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FIG. 2: (Color online) (a) The phase synchronization order 
parameter m and (b) the frequency entrainment order pa- 
rameter A/ are shown. The dimensionless length difference 
Al = \li — hi is used for the horizontal axis with h — 1, 
and each point is obtained from the random initial condition 
with Oi = and di € ±[0.05,0.1]. Clearly observed is the 
existences of two distinct dynamic phases, the antiphase syn- 
chronized state and the desynchronized state: The former is 
characterized by m ~ 1 and A/ ~ 0, while the latter by 
m significantly less than unity and A/ far from zero. It is 
to be noted that there is a broad range of bistability where 
whether the system settles down to the synchronized state 
or not depends on actual values of li and h- Change of ini- 
tial conditions is found to alter neither the dynamic phase 
transition point nor the observed bistability significantly. 



0,r = 0.8, and |/i| = I/2I = 0.35. Wc have tested the 
reduced mass fi — 0.02, 0.025, and 0.03, only to find in- 
significant differences, and the results presented in this 
work are for /i = 0.025. We use the 4th-order Runge- 
Kutta algorithm to integrate equations of motion with 
the discrete time step size At = 0.01. 

In order to measure the degree of synchrony, we define 
the order parameter m for the phase synchronization as 



m = (cos {A(j) — 7r)) , 



(3) 



where (• ■ ■ ) is the time average taken after achieving the 
steady state and A(j) = \4>i — 4'2\ with the phase ipi de- 
termined from 9i (X sin (jji. For this, we neglect the result 



generated in the initial duration of Tjni = 10^, and then 
take the average during the time T — W^. Our defi- 
nition of m gives us m = 1 if the two pendulums keep 
the antiphase synchronization, while to = —1 is obtained 
for the complete inphase synchronization. The smaller 
777, 1 is, the worse the synchronization occurs. In addition 
to the phase synchronization of oscillators, a related but 
distinct phenomenon is the frequency entrainment [ISJ. 
We also gauge the degree of the frequency entrainment 
simply by measuring the frequency difference defined by 
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where Ni is the number of oscillations of the pendulum i 
during the time T after the initial transient period Tini. 
In the numerical computation, we change I2 from 1 to 
1.006 with the interval 1.2 x lO""^, for fixed h = 1 [ll|. 
For each pair of {hjh) prepared in this way, we assign 
the quenched random number 9i £ ±[0.05,0.1] at the 
time t = and also x = i — ^1 = ^2 = is used as the 
initial condition. We then measure the two key quantities 
777, and A/ at given values of the length difference 



Al=\h-h 



(5) 



When Al = 0, the two pendulums have identical natural 
frequency and we expect them to show perfect antiphase 
synchronization Q to give us 771 = 1 and A/ = 0. As 
the length difference becomes larger, it is expected that 
beyond some value of Al the system should stop showing 
synchronization resulting in 777 < 1 and A/ 7^ 0. 

Figure [5] summarizes our main results: As is expected, 
one can clearly see the existence of dynamic phase tran- 
sition at Ale ~ 0.002 which splits the antiphase synchro- 
nized state (777 « 1 and A/ w 0) and the desynchronized 
state (777 < 1 and A/ > 0) as the length difference A^ is 
increased. The measured two order parameters are also 
shown not to have intermediate values around the dy- 
namic phase transition, indicating the discontinuous na- 
ture of the transition. Another interesting observation is 
the existence of the broad range of bistability: Whether 
or not the system approaches the synchronized state is 
not uniquely determined by the length difference Al only. 
Although not shown here, we also observe that the use of 
the larger reduced mass fx increases Ale, which has also 
been reported in Ref. 0. 

We next investigate in Fig. [3] the features of the an- 
tiphase synchronized state [(a) and (b) for A^ = 9.6 x 
10~^] and the desynchronized state [(c) and (d) for 
Al = 39.6 X 10-4] for the reduced mass n = 0.025. In 
Fig. [3)^ a) for the antiphase synchronized state, one finds 
that there is a manifest anticorrelation between 61 and 
02 due to the phase difference A<j) — it after some initial 
transient period. The approach toward the antiphase 
synchronization is clearly displayed in Fig. [31[b): Again 
after transients, Acf) approaches the odd multiple of tt. 
In contrast, the desynchronized state shown in Fig.[3][c) 
and (d) exhibit very different behaviors: The anticorre- 
lation between 6*1 and 62 becomes much weaker [see (c)] 
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FIG. 3: (Color online) The representative details of pendu- 
lum motions are depicted for the antiphase synchronized state 
[Ar= 9.6 X 10~ , (a) and (b)] and for the desynchronized state 
[Ar= 39.6 X 10"'', (c) and (d)] for ^ = 0.025. In (a) and (c), 
the time evolution of angular coordinates is exhibited in the 
plane of {61,62), green for t < 5000 and red for t > 5000, and 
in (b) and (d) the phase difference A(j}{t) is shown as a func- 
tion of the time t. In the antiphase synchronized state, (a) di 
and 62 eventually align along the line with the negative slope 
in the plane {61,62) since (b) A</> -^ n. The desynchronized 
state is characterized by (c) the scattered points in the {61, 62) 
plane, and (d) the indefinite increase of A(j> as t is increased 
due to the mismatch of the frequencies (A/ 7^ 0). 



and the phase difference A(j) increases indefinitely in time 
[see (d)], due to the nonzero frequency difference A/ ^ 
[see Fig. Hb)]. 

We finally examine that there is a hysteresis ef- 
fect around the dynamic phase transition between the 
antiphase synchronized phase and the desynchronized 
phase. For this, we fix the length of the first pendu- 
lum to Zi = 1 and change adiabatically l2- For the 
given value of I2, we integrate equations of motion for 
a sufficiently long time to achieve the steady state and 
change I2 by 6.0 x 10~^. Note that since the hysteresis 
is of interest now, the system is not reinitialized after 
the control parameter I2 is changed. As clearly shown in 
Fig. 131 our two pendulum system manifests the hystere- 
sis behavior around the dynamic phase transition: The 
phase synchronization order parameter 777 follows differ- 
ent curves when the length of the second pendulum is 
increased (the red curve in Fig. U) and decreased (the 
green one in Fig. S]). All the observations, i.e., the dis- 
continuous nature of the dynamic phase transition, the 
broad range of bistability, and the strong hysteresis ef- 
fect, strongly suggest that a subcritical Hopf bifurcation 
takes place there if stated in the language of nonlinear 
dynamics [iq . 

In summary, we have numerically studied the synchro- 
nization of the two nonidentical pendulum motions, piv- 
oting on a common movable frame. Within the limitation 
of our setup of numerical experiments, it has been clearly 
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FIG. 4: (Color online) The phase synchronization order pa- 
rameter m is measured as the length h of the second pendu- 
lum is increased (the full red curve) or decreased (the dotted 
green curve) adiabatically, while the length of the first pen- 
dulum is fixed to Zi = 1, and I2 = 1-f £2. The hysteresis effect 
is clearly seen around the dynamic phase transition splitting 
the antiphase synchronized (the upper branch of the curves) 
state and the desynchronized state (the lower branch). 



shown that the system exhibits the discontinuous dy- 
namic phase transition from the antiphase synchronized 
state to the desynchronized state as the length difference 
is increased. We have also shown that the discontinuous 
nature of the transition is reflected to the broad range 
of bistability and also to the existence of the hysteresis 
effect. We believe that our model reproduces the odd 
sympathy C. Huygens observed in the seventeenth cen- 
tury. We finally remark that it is well-known that the 
dynamic property of nonlinear system is significantly af- 
fected by such mechanical details if involved in the non- 
linearity [13]. In the present system, it is the escape- 
ment mechanisms [7|, |8| or the structural design of the 
system [l3, lll|, for example. Therefore, it is still neces- 
sary to examine the various nonlinearity setups for the 
sound understanding of the synchronization phenomenon 
in general as a dynamic phase in statistical physics. 
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